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Abstract. - We show that in a two dimensional electron gas (2DEG) the interplay between the 
Rashba spin-orbit coupling and a potential barrier {eg., the sample edge) generates interesting 
spin effects. In the presence of an external charge current, a spin accumulation is built up 
near the barrier while in equilibrium there is an inhomogeneous spin current which is localized 
at the barrier and flows parallel to it. When an in-plane magnetic field perpendicular to the 
barrier is applied, the system also develops an inhomogeneous charge current density. These 
effects originate purely from the structure of the eigenstates near the boundary. 



Spintronics has recently emerged as an active field of research due to its potential impact 
on both the design of future electronic devices and quantimr computing [1]. The ultimate 
goal of spintronics is the manipulation and coherent control of the electronic spin degrees of 
freedom. This requires the ability to generate, inject and control spin-polarized currents in 
electronic devices — an example of that is the spin-transistor proposal by Datta and Das [2]. 
In recent years, the spin-orbit coupling has been recognized as an efhcient tool to manipulate 
the electron's spin leading to a substantial amount of work devoted to study its effect on 
the transport properties of nanostructures and two dimensional electron gases (2DEG) [3-13]. 
An important step along this line was the recently reported observation of the predicted spin 
Hall effect (SHE) [14, 15]: when a transport current flows through the system the opposite 
spins accumulate at the lateral edges of the sample. It has been argued that the SHE can be 
an extrinsic effect when it relies on impurity scattering [16, 17] or an intrinsic effect when it 
originates from the spin-orbit coupling in a 2DEG [18, 19]. 

In this Letter we show that in a ballistic 2DEG with Rashba spin-orbit coupling [20,21] the 
presence of a potential barrier in one direction leads to a current induced spin accumulation 
at the sample edge. The effect is reminiscent of the intrinsic SHE although in our case it is 
based purely on geometric properties of the clean systems. In addition, we show that there 
is an equilibrium spin current [22] — ie. with no external fields applied — that flows parallel 
to the barrier, and an inhomogeneous charge current density when an in-plane magnetic field 
is applied. These intrinsic effects are relevant for a better understanding of spin-sensitive 
transport experiments and numerical results in small systems [23-26] . 
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Fig. 1 - Schematics of the band structure. The dots show the spin orientation in the two bands. 
The others panels show all possible configurations (for e > 0) of the incident plane wave and the 
corresponding spin orientation together with the two Fermi surfaces: b) incident plane wave in the 
'+' band; c) incident plane wave in the ' — ' band with \ky\ < k+, and d) with \ky\ > k+. In the latter 
case, the reflected component on the '+' band is an evanescent mode whose spin points out of the 
plane. 



We consider a 2DEG with Rashba spin-orbit coupling. The Hamiltonian is then given by 

H = + V{r) + ^ {pyo,, - p^ay) , (1) 

where V{r) is the confining potential in the plane of the 2DEG, m* is the effective mass and a 
is the Rashba coupling constant [6,27]. For a system with translational invariance {V{r) = 0), 
the eigenfunctions of Hamiltonian ^ are given by 

with e"^ = {ky — \kx)/k and A the system's area. Notice that if we write k = A;(cos(/3, sin(/?), 
then (j) = (f — n 12. Consequently, these modes have the spin in the plane of the 2DEG pointing 
in a direction perpendicular to the wavevector k as schematically shown in Fig. ^ The spin 
degeneracy is lifted by the spin-orbit coupling leading to two non-degenerated bands with the 
following dispersion relation: £±{k) — h^k^ /2m* ± ak. 

Let us now consider a system with a sharp edge in the a;-direction defined by the potential 
V{r) = VoQ{x) with Q{x) the step- function. For Vq :s> Ep, where Ep is the Fermi energy, 
the boundary condition imposes the eigenfunctions with e± < Ep to be zero for x > (the 
exponential tail inside the barrier is negligible) . For a; < they can be written as superposition 
of incident and reflected waves and for positive energy e we have: 




where s = ±1 is the band index of the incident plane wave with wavevector [k^, ky), C is a 
normalization constant, {—kx,ky) and {—k'^,ky) — k^s{cosLp' ,sanp') are the wavevectors of 
the reflected plane waves, fcs = [(2fi^e/m* -|- a^)^ — sot]m* /h^ and 

^ ^ cos([v? - if'] 12) icos^ 



sin([(p -I- (/?'] /2) ' sin([(^-|-(^']/2) ■ 

For <p 0, the reflected wave has a non-zero amplitude in both bands as illustrated in FigQ] 
(this is a simple way to create polarized beams out of an unpolarized one [28,29]). When the 
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incident plane wave is on the s — ~1 band, the above solution remains valid only if | fcj, | < 
For \ky\ > /c-)_ there are no propagating modes in the "+" band and evanescent modes localized 
at the boundary appear (Fig. Id). In that case k'^ is replaced by an imaginary wavevector and 
the corresponding spinor acquires an explicit projection on the z axis: 

*-<o ^ ^ (td ---- (:-,::) c)} . w 

with K ^ {kl~ fc2 )2 and a^/b^ = (ky + K)/{ky - K). Since the incident and the reflected 
plane waves have different spin projections, their interference leads to a non-zero spin density 
on the z-axis, {az)k,s = s('")<5'z^fe,s(''') with some spin accumulation at the sample edge 
(x < 0) . The sign of the spin accumulation depends on the sign of the conserved component ky 
of the wavevector. Consequently, in the ground state where states with positive and negative 
values of ky are occupied, the contribution to the local magnetization of each band is zero. 
However, if there were a preferential motion of the carriers in the y-direction — imposed by 
an external current, for instance — there would be a net magnetization at the surface. To 
illustrate this effect, we consider an ideal ballistic system where the voltage drops at the 
contacts and the electric field is zero at the interior of the sample. Then, in the energy 
interval [E-p + Ep — ey/2], only states with ky > are occupied within the sample [30] 

and we obtain (for Ep > eV/2): 

{a,{x))^ {^z)k,sFkAEF.eV) (5) 

k,s;ky >0 

with Fk,siEp,eV) = e{Ep + eV/2 - £,(fc)) - Q{Ep - eV/2 - e,(fc)) [31]. For small eV only 
states at Ep contribute to the sum and we may separate the contribution of the states with 
evanescent waves (Eq. (@J) from the rest (Eq. Q). For any Ep , the former contribution can 
be scaled onto the Ep = result (shown in Fig. 2a) when plotted as a function of (fc^ — fc^) 2a;. 
From this scaling law, it can be shown that the integrated spin accumulation due to this 
contribution decreases with Ep as {Ep + cP'm* I2Y? )~4. As we show below, states with 
reflected waves in the two bands (Eq. (3)) tend to cancel the effect, leading to a faster decay 
of the spin accumulation as a function of Ep. While the eigenfunctions presented above give 
a clear physical picture of the edge effect, an analytical calculation is not simple since: a) the 
two modes of Eq. (|3J) are not orthogonal to each other; b) a proper counting of modes requires 
to work with a system of finite width La,. In what follows we present results obtained by a 
numerical integration of the Schrodinger equation using finite differences. One advantage of 
the method is that it is not restricted to a square well confining potential — our results hold 
in the case of a sharp parabolic potential — and is equivalent to use a tight- binding version of 
Hamiltonian Q [32], 

ncr <n,m>(T 

(4T^(»+5)i + '^rUC(n+S)T) - {(i^C<n^x)l " ci^C(„+2)-f) } + h.C. (6) 

n 

here cjj^ creates an electron at site n with spin az — <J and energy = At, t — /i^/2m*ag, 
Oq = 5nm is the lattice parameter, and A = a/2aQ. The summation is carried out on a square 
lattice where the coordinate of site n is r„ = n^x + Uyy with x and y the unit lattice vectors. 
The hopping matrix element tnm — t connects nearest neighbors. All quantities presented 



4 



EUROPHYSICS LETTERS 






~ 6 






" 4 


A 




- 2 






_ 






-10 -8 -6^ -2 2 AAflAir 




la) 




I.I. 



> 4 
(D 

E 

c, 2 
E 


-2 



. 6x1 0" 
"4x10' 
-2x10' 



r 


' 1 ' 


' 1 ' 












-10 


-6^-2 2 Jjk 


1 . 1 . 1 . 1 . 



X [|.im] 

Fig. 2 - Current-induced spin accumulation, x^i^, Ef), close to the potential barrier (edge) for Ef — 
(a) and -Ef = 5meV (b). The insets show the convolution with a Gaussian with an rms of 0.5 fim. 



below are obtained from the one particle propagators (for details see Ref. [11]). We first 
present results for systems with a large width and then we discuss the effect of the lateral 
confinement in a quantum wire of reduced Lj. . 

To analyze the spin accumulation we linearize Eq. Q and define Xzix) — {(yz{x))/eV. 
The total spin accumulation for E-p = is shown in Fig. 2a. The inset shows the convolution 
with a Gaussian with an rms. of 0.5 /im. For E-p > (/c+ ^ 0), the contribution from the 
states with \ky\ < fc+ partially cancels the effect as shown in Fig. 2b. Thus, the maximum 
spin accumulation is obtained at Ep = and as the carrier density increases it decreases and 
tends to zero. Xz{x) is strongly dependent on Ep so even for moderate values of eV, non-linear 
effects are observed and Eq. ^ should be used. 

In this ballistic regime, the charge current also generates a local magnetization {ax{x)) in 
the x-direction {i.e. in the plane of the 2DEG and perpendicular to the sample edge) [24] 
while by symmetry {cTy{x)) = [4] — the non-zero value obtained for finite systems [24] is 
a consequence of the reflections at the sample-lead interfaces and depends strongly on the 
characteristics of such interface [33] . It is important to emphasize that the spin accumulation 
at the edge is a purely geometric effect. It is not directly connected to the one discussed in 
Ref. [18]. The latter is due to an electric field induced spin current, which is absent in the 
ballistic regime studied here. 

Another consequence of the symmetry of {az)k,s with respect to ky is the presence of spin 
currents along the edges of the sample. In equilibrium, we define the spin current density as: 

h 

occ occ 

here iiy — i/h[H^y\ — py/m* -I- aax/h is the y-component of the velocity operator and the 
summation is restricted to all occupied states with quantum numbers ky and v. Here again 
the evanescent modes play a central role and Jy'(x) is non-zero in equilibrium. Figs. 3a 
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Fig. 3 - Total spin (a,b) and charge (c,d) current density in the vicinity of the edge for different values 
of Ef. Notice the currents are localized near the edge. Here qubBx ~ 0.25meV and a = 5meVnm. 

and 3b show the total spin current density Jy'{x) for different values of Ep. The total spin 
current oscillates with a characteristic wavevector fcF,+ + kp.- and decays toward the bulk of 
the sample with a characteristic length (/cf,- ~ fcF,+)^^ — /2m*a k, lOOnm that is energy 
independent. This longitudinal spin current, integrated in half of the sample — from —Lxl2 
to — is non-zero: spin flows in opposite directions on opposite sides of the sample. It is 
worth pointing out that this current does not violates time reversal symmetry [22]. 

The counterpart of the current induced magnetization at the sample edge, is the gener- 
ation of edge currents by a magnetic field. To avoid the diamagnetic coupling and Landau 
quantization, we now analyze the response of the system to an external in-plane magnetic 
field. In the continuous model, the propagating modes are still given by Eq. Q but the 
spinor's angle is such that e"^ = z/\z\ with z = a{ky — ik^) + \g^B{—Bx + i-By) and g the 
g-factor. The presence of the magnetic field modifies the Fermi surface which is now given by 
Ey = h'^k'^/2m* ± \z\. Notice it is not longer a circle. It is then clear that the magnetic field 
introduces an asymmetry between ky and —ky when applied in the x-direction. In bulk, the 
charge current induced by the field is zero due to a cancellation between the contributions 
of the two bands. However, at the edge the evanescent modes lead to the appearance of a 
non-zero charge current density, 

Jy{x,e) = {vy)ky,u5{e- Ek^^y) (8) 

The total current density, Jy(x), is obtained by integration of this quantity up to E-p. Figs. 3c 
and 3d shows Jy{x) for E-p = and 5meV. As in the case of the spin current, the charge cur- 
rent is localized at the surface and presents modulations with two characteristic wavelengths, 
Af/2 = 2Tr/{kF,+ + kp,-) and Xa = 27r/(fcF,- — kF,+)- 

So far we have presented results for wide (large L^) systems. We have shown that out of 
equilibrium, when a ballistic transport current flows in the system, there is spin accumulation 
at the edges. The spin and charge current densities discussed above are properties of the 
ground state and depend on a delicate balance of the contribution of states with and without 
evanescent waves. At a given energy, the number of states of each type can be changed by 
engineering the sample geometry and thus altering the partial cancellation. In particular this 
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Fig. 4 - Current-induced spin accumulation (a, b) in the absence of a magnetic field and charge (c, d) 
and spin (e,f) current density in the presence of an in-plane magnetic field for a quantum wire. The 
panels on the left (a, c and d) correspond to the case of a non-resonant energy, e = 4.9meV, while 
the panels on the right (b, d and e) correspond to resonant cases, e = 4.98meV, e — 4.95meV and 
5.04meV respectively. We used: q^bBx ~ 0.25meV, L = 500nm and a = SmeVnm. Panel (g) shows 
the density of states in the presence of the magnetic field. Notice the singularities are split. 

can be done in a narrow (small L^) quantum wire. The lateral confinement of the 2DEG in 
the x-direction leads to the quantization of the transverse modes and the density of states 
presents quasi-one dimensional van Hove singularities (see Fig. 0^) . When E-p is tuned to be 
close to a singularity, there are many states with ky ~ — thus with no evanescent waves — that 
define the behavior of the system. 

The results for the current-induced spin accumulation, obtained using the tight-binding 
approach, are shown in Figs. ^ and^ in the absence of a magnetic field. Since the wavefunc- 

tions of the states with fcj^ — have a simple form ^'j- (a;) = Lx^ e'}cp{^mx/Xa)sm{mrx/Lx) (^'), 
where n is the channel index, one can readily show that, for ky > 0, (az) oc sin^ (nnx/Lx) 
sin(27r(a; — Lx/2)/ Xa)) in excellent agreement with the numerical result of Fig. 4b (n = 15 
for Ep « 5meV). In this case, the sign of the magnetization at the edges also depends on the 
ratio Lx/Xa- For the parameters used in Fig. 4b the spin density is positive in the interval 
[0, £3,72]. Figs. ^ and^Ji shows the spin current profile. The resonant case presents well 
defined signs: there is a net spin current flowing in one direction on the right and in the 
opposite on the left. 

The charge current density in the presence of an in- plane field is shown in Figs. and 
0f. For \Bx\ ^ 2akp/giiB a similar calculation as above gives {vy)± cx -iisir? {nirx/ Lx) 
cos{2tt{x — Lx/2)/Xa))- Though this expression does not fit the numerical result exactly 
(Fig0f), it contains the correct spatial modulations. In particular, it shows that the charge 
current density changes its sign in a length scale given by Aq/2. It is also straightforward to 
check that the van Hove singularity is split by the magnetic field — the energy splitting being 
A « gusBxl siii{nLx/Xa)\Xa/LxiT for n > Lx/Xa (see Fig. 4g). 
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In summary, we showed that the presence of a potential barrier in systems with spin-orbit 
coupUng leads to the appearance of current-induced magnetization at the sample edges and 
of equilibrium spin and charge current densities. These effects are important for a correct 
interpretation of numerical results in finite systems, where geometry plays a crucial role. On 
the experimental side, low density ballistic 2DEGs or quantum wires are required in order to 
observe the spin accumulation described here. It is worth mentioning, that these conditions 
are not met by the recent experiment of Ref. [34], where disorder is very large. 

^ ^ ^ 
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